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NEUTRON PHYSICS PARAMETERS OF SOME WATER-MODERATED 
URANIUM LATTICES 

LG.AndreeVy LII.Ganev, S.I.Garanijij EjkDvoinishnikov^ A.I.Dokutchaev, F.S. Drozdov, A. F. Zaitsev, 
V.P. Katkov, L,V .Komissarov, J u.M.Knisznikov, JuJ a. Kravtchenko, N .A.LasiLkov, Ju.V. Nikolski, 
A.N.Protzenko, G.A.Stoljarov, V .M.Talizin, l.A.Filatjev, Ja.V.Shevelev. 

The results are reported of three independent groups of experiments on light water lattices 
with rod elements. Some of the lattices considered (those incorporating no absorbers) are 
rather similar to ones described in the publication ''Light water lattices", edited by IAE]A 
in 1962, with some differences in enrichment and rod diameters. 

In section A the experiments arc (b^scribed on 5% enriched uranium dioxide systems with 
fuel rod diameters 0.45 cm. Furthermore, r>ift>ry lattices are considered in this section, includ- 
ing uranium and absorber rods, with optical transparency in the range of 0.75 to 40. The 
results are reported concerning the critical masses of the lattices and some relations are 
presented which has been used in calculations. 

In Section B the multiplication parameters of the natural uranium dioxide lattices are 
studied by the method of absorption ratio measurement for core components, the neutron 
leakage correction being taken into account. 

In section C performances of the lattices with magnesium — 2 % enriched uranium alloy fuel 
rods are studied by the pulsed neutron source technique. 

This study is performed in Kurchatov Institute of Atomic Energy. The authors of the 
separate parts of this report are; 

Section A - l.H.Ganev, E.A.Dvoinishnikov, A.I.Dokutchaev, A.F.Zaitzev, i\ .A.Lasukov, 
A.N.Protzenko. 

Section B - S.l. Garanin. V.P. Katkov, Ju,N.Knisznikov, L.V. Komissarov, Ju .Ja. Kravtchenko, 
Ju.V. Nikolski, G.A.Stoljarov, l.A.Filatjev. 

Section C — L.G. Andreev, F.S. Drozdov, V. M.Talizin, Ja.V.Shevelev. 

A. CRITICAL EXPERIMENTS WITH LATTICES CONTAINING ABSORBER RODS 

One-group method was used in material buckling calculations for binary lattices. Epithermal 
processes were taken into account by additional factor, introduced into the infinite lattice 
multiplication constant, according to ref. 1. A considerable reduction in computer time results, 
which makes possible the precise calculation of the spatial non-uniformity and the time 
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dependence of the lattice properties due to burnout. The disadvantage factors were calculated 
in the diffusion approximation, the effective boundary conditions being imposed on the rod 
surface, A formula was obtained on the basis of consideration of theeffective boundary con- 
ditions for absorbing cylinder in an infinite non-absorbing medium, realized in ref. 2. This 
formula does not take into account the scattering in the inner region of the rod and allows 
the analytical integration when averaging y over neutron spectra 


y- nod - e foe- 


q 2. 


(A-1) 


where 


fi « y — ; a = d 2^ 


/8 =dv^ (294°K) 


0.0253 

E 


The values of the coefficients h^, h 2 and h 3 for various i are given below. 

1^234 5 6 8 9 

1>1 0.7057 0.3 0.5526 0.71 0.1974 5.52 0.8281 0.2367 

h2 0.085 0.243 0.5908 0.32 0.392 0.1765 ’ -1.2934 -0.3069 

h3 2.58 1.337 1.748 1.46 1.2253 0.6044 2.77 2.77 

The quantities y and g^jj were averaged over the Maxwellian spectrum when imposing the 
boundary conditions and calculating the self-shielding factor for strong absorber in the thermal 
region . 

The self-shielding factor was calculated from the neutron balance condition, taking into 
account neutrons, incoming on the rod surface and absorbing in it, 

g L (E) = 4^ ^ 1 

0.75d2^(E)7 (A-2) 

On averaging over the Maxwellian spectrum, we obtain 

[^ ( Va e -> - ifi Hi ^ [1 1— 3/ J , ^ , 

^ ^ a+V/sV (Uf4//32)2 


(l + f6//82) 


__J [e-x(l^ 

■v/>r'0.75 


[e-x(l^) _ g-a /2(_l+2£x 

\fi'x 


foA 

ai=l 

^0 


/3o = i9(E = Eth). 


(A-3) 
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It was found convenient, in calculating the average self-shielding factor gg for epythermal 
region, to use the foliOwing expression for the quantity y; 




(A-4) 


The coefficients fg and fp are to be found from the equation for fj with the above values 
of hj, h 2 and h 3 ; the coefficients fy and fjQ are determined by Eq.* (A~5). 


^3 +f5 


2(0.75 -fyfe) 


On averaging over the slowing-down neutron spectrum we obtain: 


J___ rr^ ir, ^ To /3 j \ ^ 


(A^5) 


oWs ) -► flO *fP33)l (A-6) 


iSs - A 


, '\/^th y Ps “■ 

'oV-fr 5 ^ 

l-S \/l+y^ 2 + fQiSq 


EXPERIMENTAL ARRANGEMENT AND RESULTS 

The core was composed of rod-type fuel elements. The average length of the fuel element 
tube, filled with sintered uranium dioxide pellets, was equal to 59, 66 cm. Steel- and zirco- 
nium-clad fuel elements with the inner diamete" of 0.45 cm.were used in the experiments. 

5.32% and 5.12% enriched uranium was contained in steel- and zirconium-clad fuel elements, 
correspondingly. The average weights of the fuel element cores were 76.51 g and 91.63 g, cor- 
respondingly. The absorber rods used in the experiments were constructed of steel tube filled 
with the B 4 C -f- AI 2 O 3 mixture on the length equal to the core height. The boron content in the 
boron carbide amounted to 65.6 weight percent. 

In the first series of the experiments “pure" assemblies were studied, containing no boron- 
filled rods. The rods were arranged in a triangular lattice, the lattice pitch being 7, 8 , 14 mm 
for steel and 8 mm for zirconium-clad fuel elements. 

Besides the "homogeneous "assemblies described criticality measurements were made on 
the assemblies with periodical core heterogeneity. 

Some typical lattice cells are illustrated in Fig.A -1 (the boron rods shown in Fig.A— 1 were 
used in the second series of the experiments). 

The assemblies were surrounded by practically infinite water reflector (~20 cm.). The 
water temperature in various assemblies varied in the range of 20 ± 3 ^C. 

The experimental and calculated values of the critical mass for various s pacings of 
"homogeneous" and “heterogeneous" assemblies are listed in Table A-1. The uncertainty 
356 
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Critical masses of the '*pure'^ assemblies 

Table A— r 


Fuel 

element 

type 


Lattice cell 
configuration 

^0 

SU02 

^11 

P5 



.,culc 


7 

’MIomogeneous** 

1.53 

100 

14.0 

3937 

15.1 


8 

'Mlomogeneous'' 

2.39 

156 

7.6 

2124 

8.1 


8 

87-1, no gap 

2.49 

163 

7.2 

2030 

- 


8 

37-7, no gap 

3.25 

213 

5.8 

1626 

- 


8 

37, single row gap 

3.52 

230 

5.6 

1589 

5.5 

Uclad 

8 

37-1, single row 
g«P 

3.62 

237 

5.5 

1538 

5.5 

Stee 

8 

37-7, single 
row gap 

4.61 

301 

4.7 

1335 

4.8 


14 ' 

’’Homogeneous'’ 

10.06 

658 

4.7 

1319 

5.1 

3 

’S "2 

8 

’’Homogeneous” 

1.64 

98 

8.8 

2138 

- 

oiS 

8 

37-1, no gap 

1.73 

104 

8.1 

1974 

- 


8 

37-7, no gap 

2.45 

147 

5.7 

1380 

- 


Critical number of fuel elements and number of boron rods in two-zone assemblies 

Table A-II 


No. of 
experi- 
ment 

Lattice 

spacing 

mm 

Central zone 

Peripheral 

zone 

cell 

Number 
of boron 
containing 
cells 

Critical 
number 
of fuel 
elements 

Lattice cell 

absorbe 

r rod 

boron 

content 

g 

cladding 

din. 

mm 

1 

7 

61-1, no gap 

7.35 

:v . ■ " 

I’MIornoge neons 

” 28 

5641 

2 

8 

37-7, single 

7.35 

5x4.6 

.>7-7, single 

77 

3461 



row gap 



row gap 



3 

8 

37-7, no gap 

0.71 

10x8 

37-7, no gap 

12 

2112 

4 

8 

37-7, no gap 

1.6 

14x12 

37-7, no gap 

7 

2119 

5 

8 

37+11-7 

25.6 

10x8 

37+11-7 

19 

5394 

6 

8 

37+11-7 

57.74 

14x12 

37+11-7 

31 

5648 

7 

8 

37+11-7 

2.78 

10x8 

37+11-7 

31 

5510 

8 

8 

37+11-7 

7.42 

14x12 

37+11-7 

54 

5636 

9 

8 

37+11-7 

1.38 

10x8 

37+11-7 

37 ^ 

5504 

10 

8 

37+11-7 

3.3 

14x12 

37+11-7 

67 

5525 

11 

8 

37+11-7 

0.71 

10x8 

37+11-7 

46 

5522 

12 

8 

37+11-7 

1.6 

14x12 

37+11-7 

93 

5527 
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of the measured critical mass values is 1.5%. The discrepancy of the experimental and calculated 
values does not exceed 10 %. 

In the second series of the experiments two-zone assemblies were studied, containing boron 
rods in the central zone. Various types of the lattice cells are illustrated in Fig.A— 1. The 
critical assemblies as well as the zones containing boron rods approximate the cylindrical 
geometry. The results of the criticality measurements of the two-zone assemblies (at T= 20 °C) 
are listed in Table A— 11 . 

The calculated and measured values of the material buckling of the absorber rod containing 
zone are plotted in Fig. A~2 against rod blackness. 

The experimental data are in good agreement with the calculated values, thus giving ground 
to use the above relatively simple technique in a rather wide range of pu / and d^p values 

P5 

for the boron rods. 

B. EXPERIMENTAL DETERMINING PARAMETERS OF 
NATURAL URANIUM DIOXIDE LATTICES 

1. PHYSICAL PARAMETERS OF THE LATTICE 


The definitions of neutron maltiplication and slowing-down parameters used in this paper 
do not imply the subdivision of neutron spectra into parts and are similar to those of ref. 3 . The 
parameters studied were determined experimentally using the technique described in the same 
paper. Furthermore, the corrections for neutron leakage were introduced in the one-group ap- 
proximation. Used in the correction calculations migration areas from generation neutrons in 
U-:^35 to fission in with taking into account cascade multiplication and absorption due 


93 Q Hz ^ 

to mam U“ resonances (M ^ ) were determined from poisoning experiments, with boric acid 

solution. It is possible to express in terms of the neutron absorption ratio for poisoned 


and non-poisoned systems, on the assumption, that poisoning does not change the neutron spectrum 
in the above resonance region: 

nf 






pois 



5 pure 




B2 - b2 • ( 5. ) /( 8 ) 

*^Dure ^Dois V r L •_ / \ — —) 


(B-1) 


^pure - opois V— t 

^ fjf pois/ f pure 

5 / ‘5 


where is buckling, and are numbers of capture events, followed by and 

O D 

fission, correspondingly. is determined similarly. 

o 

The data obtained from poisoning experiments make it possible to determine and to 

evaluate the form of the critical equation 

558 


— 5 — 
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b2(Cb) = -a [1 - e 


(Co - Cb)1 




In = (Co - CB)m , (B- 3 ) 

where Cg is the boron concentration. 

The coefficients A, Cq, 1 and m are to be determined from the experimental data. Taking 
into account, that 

1^2 _ ^ m ,13.4, 

b2^o ^ B2 = 0 


the critical equation can be rewritten as: 


,1 , . K. 

n 


Here ^ 


2. EXPERIMENTAL EQUIPMENT 

Two research reactors were used in the experiments. The first one was designed to measure 

the following parameters: g, (A. 0 . veff, P and ^ . It was a critical assembly, compris- 

5 

ing the central ’’seed” zone and the lattices investigated around it. The triangular lattices 

under study were assembled from fuel elements with spacings of 13, 15, 17 nnd 19 mm. I he 

elements used were fabricated from aluminium tubes Hlled with natural uranium dioxide 

pellets of 8.0 g/cm^ density to a height of 1200 mm. The outer tube diameter was 11.3 mm, 

the wall thickness 0.7 mm, the pellet diameter 9.4 mm. ^ 

The other research reactor was used to measure (M^) , ^nd to evaluate the form 

o 

of the critical equation for 13 mm spacing lattice. It had a form of a hollow cylinder with a 
hexagonal tank in the central region containing the lattice under study. The arrangement 
described made it possible to poison the lattice by the boric acid solution and to conduct the 
experiments on a dry lattice (without water in the tank). 


3. RESULTS 

The results of the measurements of the lattice parameters are presented in Table B-1 

and in Fig. B-1 (See Table B-Il for designations). 

The reactor of hollow cylindrical geometry has been used for the experiments with 13 mm 
' . f 9 

spacing lattice, both dry and poisoned with boric acid. On the assumption of Mg value 

being in the region from 0.6“'0.9 for dry lattice (for boron-free lattice imbedded in water the 
above ratio was found to be 0.4), the following data have been obtained for dry system: 
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nf / r =0.60+0.08; = 0.22 -0.25; /« = 1.14 -1.25 

8 

0 

(Mf )2 = 0.6 - 0.9; P = 0.8+0.3; = 0.48+0.02 

The above values are to be regarded as preliminaiy estimates and are subject to more 
accurate definition. 

The results for 13 mm spacing boron-poisoned lattice are listed in Table IB-II. 

As an analysis of the experimental data has shown, the B^(Cg) and K^^CCg) variations 
can well be described by Eqs. (B-2) and (B-3) with the following values of the coefficients: 

A=130-10-4‘; Co = -0.5; 1= 0.294; m * 0,296. 

From the data obtained one can find = 77.0±3 cm^; n = 1.0 ± 0.01. 

As it will be seen from the above results, the lattice under study can be described by 
the one-group critical ecquation within the limits of the experimental errors. 

4. COMPARISON OF EXPERIMENTAL RESULTS WITH CALCULATIONS 

Commonly used scheme was applied for the calculations of the lattice neutron multiplica- 
tion parameters; however, the values presented have the meaning adopted in this work. 

As the transparencies of the lattices investigated for fission neutrons are very high, it is 
possible to calculate ^ from the formula for homogeneous mixture making use of the published 
data on effective cross sections (See ref. 4). 

The values of g thus obtained are in satisfactory agreement with the experimental data 
(See Fig. B-1), but for 13 mm spacing lattice, where a substantial discrepancy is observed. 

When evaluating 0, equation (53.18) of ref. 2 has been used for the calculation of resonance 
integral. 

According to the results of the analysis of available experimental data on resonance capture 
in uranium dioxide blocks, the numerical values of the coefficients of the volume and surface 
terms were assumed to be 10.4 and 31 .25, correspondingly (the Vy part of the cross-section 
being neglected). 

The calculated values of 0 are in good agreement with the experimental data for all 
spacings. 

The cantities and 0 have been calculated v>^ith the use of Wilkins equation. The 
Wilkins distribution was assumed for neutron energies up to 0.5 ev, while at higher energies 
the spectrum vvas supposed to be Permian slope. The effective cross sections were obtained 
by averaging over the above spectrum, the screening effects, dependent on lattice geometry 
and neutron energy being taken into account. 

'rhe procedure described yielded the values of which are not quite in satisfactory 

agreement with the experiment, the better agreement being observed for d. 

The calculated values of neutron multiplication parameters are shown in Fig. B-1. 

7 r c 

J w u 
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Measured values of lattice parameters 


Table I3-I 


^ Lattice spac- 

mm 

Parameter 

13 

15 

17 

19 

' "s ‘ 

0.13310.008 

0.07510.005 

0.07610.004 

0.051l0.0a4 

0.05210.003 

0.044+0.003 

0.042+0.003 

ns 

1.22510.010 

1.20910.009 

1.20210.010 

1.198+0.010 


1.61110.050 

1.23810.004 

1.10010.040 

0.983+0.035 

n„/n^ 

0.22210.009 

0.47410.012 

0.72810.018 

1.010+0.025 


nAl^lg 

/ l U, 


(‘4' 4), • 

Pg - 9 5 lattice 

^ ' ^^Maxwell 

“ll20/% 


0.067+0.0at 

0.70210.015 

0.6110.03 

0.1110.01 

0.10510,003 

0.08010.003 

0,87610.025 

0.1 0910.008 
0.70910.020 
0.87010.005 
1.28210.020 
1.3210.04 

1,2910.02 

1.06310.030 


0.06710.003 

0.68710.013 

0.4710.02 

0.05810.003 

0.04 

0.14010.001 

0.87810.020 

0.063l0.0a4 

0.81510.015 

0.77810.006 

1.30310.U::. 

1.02210.025 

1.1710.02 

1.1210.03 


0.06710.003 

0.68110,013 

0.3910.02 

0.01010.003 

0.03 

0.10010.002 

0.83610.015 

0,04510.002 

0.86510.015 

0.70310.009 

1.31210.017 

0.91410.020 

1.1210.01 

1.1110.03 


0.06610.003 

0.67810.011 

0.34010.025 

0.02810.003 

0.02 

0.08110.001 

0.78810,017 

0.03810.002 

0.90610.018 

0.65810.010 

1.31710.020 

0.81810.020 

1.0710,01 

1.1010.03 


* The values of in the first line of the table were obtained in the experiments 

with depleted and natural uranium as indicators, the second line contains the values, obtained 

I f I* I* 

from cadmium ratio measurements on natural uranium and ( — - = ^ 


235 I “e - “5 


All the data listed in the table are obtained on averaging the results of two experimental 
groups. 35 g 
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Poisoned lattice parameters (13 mm spacing) 


Table 13— li 


'\I'3»03 

concentration 
mg /ml 

Para- 

meter 

0.0 

3.96 

6.35 

10.80 

14.70 


0.133±0.008 

- 

- 

0.14910.010 



1.225+0.010 

- 

1.230+0.005 

1.235+0.005 

- 

"8'"5 

1.611+0.050 

- 

1.9110.03 

2.15+0.03 

-* 


0.067+0.004 

- 

0.066+0.005 

0.06810.005 

- 

n„/n> 

0.222±0.008 


0.22610.009 

0.23110.009 

- 

n,/nj 

0 

- 

0.43+0.01 

0.7310.01 

- 

n 

CO 

0.61±0.03 

- 

0.65+0.02 

0.6810.02 

- 

® V ^ 

O.lliO.Ol 

- 

0.1310.01 

0.1710.01 

- 

J 

5 

0.105+0.003 

- 

0.09+0.01 

0.1010.01 

- 

J 

9 

0.08010.003 

- 

- 

- 

- 


0.875+0.025 

- 

0.70710.010 

0.62310.010 

- 

^-1 

0.10910.008 

- 

- 

0.10110.008 

- 

0 

0.709+0.020 


0.70010.008 

0.692+0.003 

- 

e 

0.87010.005 


0.72710.010 

0.65310.010 

- 

''eii 

1.28210.020 

- 

1.2610.01 

1.2:10,01 

- 

P 

1.3210.040 

- 

1.55+0.02 

1.73.;0.03 

- 

(M^) (cm^) 

8 

33+5 

- 

- 

- 

- 

8 

8015 

- 

- 

- 

- 

b 2 • lo4(cm-2) 

16.8+0.06 

30.710.05 

37.510.08 

47.410.16 

53.3+0.14 


Here Hj/rij is neutron absorption ratio for i and j— th components of the lattice; 
a- - is the share of i— th component in epycadmium neutron absorption. 

Indices refer to Al, 11 and !3 in the core, indices 

f, r, R designate fission, radiative capture and resonance absorption on correspondingly; 


Cy and are cadmium ratioes for and 

256 
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C. UTTICE PARAMETER STUDY DY PULSED NEUTRON SOURCE TECHNIQUE 

A common method of studying the properties of neutron multiplication systems implies 
the measurement of stationary neutron density distributions obtained in critical experiments 
or on subcritical systems due to constant neutron source. In many instances, however, more 
precise results can be obtained in a simplier way by studying non-stationary distributions, 
produced by initial neutron pulse. The quantity measured directly in this circumstances is the 
time constant of neutron density decay at a given point in the system. Fven more effective 
are the experiments, employing both the methods simultaneously. 

Both in planning the experiments with a pulsed neutron source and treating the obtained 
results the approximate relations are of use given below, describing space-time distribution 
of neutron flux. The amplitude of the spatial harmonic associated with buckling in 

an interval between two pulses varies according to the relation: 


(t) ~ exp (Snot) + — -- 


Pn ^no ^ 


(C-l) 


where is reactivity (expressed in terms of ^gff) corresponding to n ~ th harmonic, i.e. 

1 -[K„T(b2)]-1 1 -7(132^,.)/ ((13-) 

Pn= rZ = * (C-2) 




The function f(B^) is the leakage escape probability i.e. Fourier transfonn of the second 
neutron generation sources distribution for point primary source in infinite media. If 
Koo — 1 > /Seff Pn negative, when n ?^0. 

The parameter 

X T f(B2^,)/T(i;. 


(1-/5,) 


is very small, since the first co-factor is very small (A = I Aj is the mean \alueof the 

i fi 

decay constant for delayed neutron sources, T is the prompt neutron life time), and the second 
co-factor in any case does not exceed unity, 3“ designates the value for a system of 
critical dimensions. We have for 


(l-Pu)/^ 


“o (1 + (i\) = 

3“ ) 

‘"^0 cr 


^ ““ ^effPn 


Pen 

(U.,). (C-1) 


The interval r between two pulses is supposed to satisfy the following inequalities 

iSnol r » 1; A r « 1 

both being possible to hold simultaneously, since | similarly to f, is very small 


(C-5) 




10 — 
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As it follows from Eq. (C-1) and inequalities (C-5), neutron flux c;i!>(t) can be regarded as 
a sum of rapidly varying part and '’background”. In turn, the background is a sum of a constant 
and approximately linearly decaying term. The to^l background variation is only a small part 
of the rapidly varying component, of the order of I'he Ixickground-to-rapid component 

ratio depends upon the degree of the system subcriticality and the pulse repetition frequency. 
When I SQ^rl is supposed to be large and constant, it follows, that the background increases 
in inverse proportion to (being lower for n > 0 harmonics). The background to the most slow 
exponent amplitude [i.e. exp (Sp^t)] ratio for critical system increases after each pulse, the 

increment being T/ _ , where 7 = 5— . When the number of pulses totals to /Seff^- 

//Seffr 1 

or at a high pulse repetition rate even before (when a total number of pulses ^eff/^ ^ 

attained), the background becomes equal to the initial value of the rapid component, and the 

system is to be shut down for a long time. With I /DqI ^1 the rapid component already can be 

separated easily from the background. The main exponent exp (S^^t) begins to dominate after a 

short time, since for nearly critical system 


Soo 1 - po f (B^) 


(C-6) 


The increase of higher harmonic contribution in strongly subcritical system makes it 
difficult to measure the quantity . When I p^l » 1 we obtain, instead of Eq. (C— 6): 


l-T(B2)/f 

l-f(B2)/r(B2„) 


(C-7) 


In the limit 0 the above ratio becomes unity. Without having resource to the special 
methods of the source-detector disposition, it is possible to determine at f (B^) > 

The measurement of the Sqq versus system dimensions (B^) dependence makes possible 
the determination of important media parameters. E.g., near-critical experiments yield a 
quantity 


^eff cr 


(C-8) 


the measurements, performed in a wide interval of B^ values offer the possibility to determine 
the forms of the functions f (3^) and K^. Extrapolation of f (B^) to = 0 yields 


Kco = 1 - i^eff 


SqoCO) 

Soo<b2) 


,2 _ WSQQ/dB^)g2_Q 


(C-9) 
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The above relations imply that the time of birth of any next generation neutron and its 
distance from the source are not correlated statistically. This assumption is well satisfied 
in water-moderated uranium systems, since the neutron migration is small (L^ « M^) in this 
case in the energy interval of significant absorption. Furthermore the dependence of /5^ff on 
system dimensions has not been taken into account. The corresponding corrections are found 
to be very small. 

The following relation can be obtained for zeroth harmonic: 




„9 . S„„(b2) 




(i+r' 4 G+ii+j) . (c-ii) 


Here F, G, II, J are corrections for thermal neutron diffusion, capture during slowing down and 
thermalization, deviation of the capture and fission Sf cross sections from Vy law, and 
blocking-effect, correspondingly. When L“«M- , we obtain 


2 9 9 

F={l-_2 )B L“ ; L~ =!J T 

r) OCrcc oo eooo 


(C-12) 


(the infinity index refers to the neutrons in thermal equilibrium with media). 


2 = F/ tSD(t) 
o 




; SD = D(t)-D 


(C-13) 


Here D(t) is diffusion coefficient averaged over neutron spectra which forms at a time t after 
neutron birth 


H = H..; 


,, ,^00 ,, ~ 5 (Sfv) dt 

epith - (- — — ; 

^ocr ° (^fv)oc 


(C-14) 


S(2f v) = (Sfvlt - (Sfv),„ . 

In particular, the epithermal contribution to [1 is determined by 


‘^epith ~ (1 - 7 ) ^eff 

^00 cr 


Jf S 


a * - «> 


F" ' ^fff. • 

^"^s Ifv 0.1 ev E 


(C-15) 


The heterogeneity correction can be evaluated approximately from the following formula: 

~ ^^00 ^00 cr^ 




'-J1 - 0) 


T|„ (1 - 0) + 'I'f [i9 + (1 + / . ) (1 - 0 )] ’ 
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where Tf, ^^re the life times of neutrons in fuel rod and moderator, are coirespond- 

ing fluxes and 0 is neutron absorption probability in rod. 

In buckling (B^) calculations it is important to take into account the extrapolation distance 
(5) dependence on the core boundary curvature. Furthermore, 5 is influenced by the non-statio- 
nary character of the process itself and the value of Sqq = a. In fact, the non-stationary 
equation with exponential decay is an equivalent of stationary equation with the effective 
capture cross section 

2.fr - 2 • (c-17) 

However, S is dependent on the capture cross section. Thus, decrease in S^ff, obviously, 
would result in increase of B, The physical explanation of this fact follows from the considera- 
tions given below. The neutrons, diffusing in the reflector core at any time, have arisen some time 
ago, when the neutron density in the core was higher than presently existing. Hence, the 
core reflux to leakage ratio at any time is higher than in stationary conditions. That is, the 
reflector efficiency during the neutron density decay is higher than in stationary conditions. 

Cadmium shielding of core is a radical method of decreasing variations of 5. Then the 
reflector efficiency is determined by fast neutron reflux only, i.e. it does not depend on the 
effective capture cross sections and, hence, on a. The dependence of 3 on the core boundary 
curvature radius was proved to be weak if the core and reflector diffusion coefficients are ap- 
proximately equal. 

Without resorting to stationary neutron distribution measurement, the values of B and the 
dependence of this quantity on core boundary curvature seem possible to be determined from 
experimental a vs.B^ curves, obtained for a wide range of variation of the factors B depends 
upon. For example, comparison of values obtained for cylindrical cores of different height 
and radius but characterized by the same value of a offers the possibility, in principle, to 
determine extrapolation distance. 

The authors have studied water-moderated uranium system, containing 2% enriched uranium- 
magnesium alloy rods arranged in triangular lattice of 60 mm spacing. Neutron source and a 
thermal neutron detector were located in the core. The zirconium tritide target (8 mm dia.) of 
the neutron generator was used as a neutron source. A scintillation counter with Li^F + ZnS 
powder, pressed on the organic glass light pipe butt-end was employed as a suitable detector 
satisfying main requirements: minimum neutron field perturbation and high sensitivity to 
neutrons combined with good neutron- y-ray discrimination. 

The measurements of thermal neutron density time dependence have been performed for 
12 core configurations differing in number of fuel channels and height, though all the assemblies 
being nearly cylindrical in shape. The thermal neutron density time decay constant a was 
measured versus buckling for assemblies of two types: 1) cadmium shield-free assemblies 
surrounded by water reflector, 2) cadmium shielded assemblies, i.e, ones separated from 
reflector by cadmium sheet. The cadmium sheet dinged to the core tightly, reproducing 
elementary cell configuration on the core boundary. Specifical positioningof the source (at 
the core center) and detector (at 2/3 of core height and 0.44 of radius produced the most 
356 
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favorable conditions for singling out the fundamental mode in neutron distribution (See 
Fig. C-1). 

The values of the axial (5{-|) and radial (5j.) extrapolation distances were determined from 
radial and axial neutron flux distributions mea'iUred with indium foils in critical assembly. The 
radial flux distributions have been measured along five radii in the central plane. 

To prove a weak dependence of the decay constant a on the deviation of the core from 
cylindrical shape check experiments have been performed on one of the assemblies with a 
small number of fuel channels. The deviations substantially exceeding those inherent in the 
parameter measurements were found to produce no significant changes in a. Hence, the use 
can be considered for buckling calculations of the following ideal cylinder formula: 


_ 2.405 


H +25|| 


(C^18) 


where A : 


is the migration area anisotropy ratio. 


1 

A comparison of the family of curves a = f(B^) obtained for core configurations of the 
same radii to those corresponding to different radii made it possible to determine anisotropy 
ratio A and dependence of on core boundary curvature radius. For given values of 5|| and 
Sro (^0 5r value for critical assembly) the above curves can be matched at definite A 
values and definite form of - f(r) only, the deviations of the curves due to these two 
factors displaying different behaviour and being sensitive to the change of A and 5^, = f(r). 

The curve matching at the values of and obtained experimentally yields A = 1.1 and 
displays no dependence of 5 j. on boundary curvature. The following good approximations were 
found by the least square method under these conditions: for assemblies surrounded by water 


= - 1954 + 706000 - 21400000 


(C-19) 


for cadmium shielded assemblies 

a = -1928 + 664000 - 1 1400000 B^ {C-20) 

illustrated in Fig. C-2. The difference between Eqs. (C- 19) and (C-20) could possibly be 
attributed to the influence of the non-stationary character of the process on extrapolation dis- 
tance in the cadmium-free assemblies. The difference decreases with the increase of core 
dimensions (the influence of the non-stationary nature of the process becomes less significant 
as a decreases). The Eqs. (C-19) and (C-20) yield nearly equal values of 

According to the three group calculations 8^ does not depend on the boundary curvature 
radius for cadmium-free and cadmium shielded assemblies. The calculated extrapolation dis- 
tance for cadmium - free assemblies coincides with the experimental value, whereas for 
cadmium shielded assemblies it lies higher. The maximum difference between measured and 
calculated buckling values, the latter including the calculated correction for non-stationarity, 
amounts to 1%. The corrections for thermal neutron diffusion, capture during thermalization 
and diffusion, deviation of the cross sections Xg andXi* from law and blocking-effect, are 
insignificant (1%) and of different signs. 
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The parameters of the described water moderated uranium system determined in the present 
work, are listed below. 

1. Critical load at 22°C (core height 80 cm) 

N^r = (204.9±0.2) fuel elements 

2. Material buckling 

b 2 =r 3.229. cm-2 
cr 

3. Buckling for prompt criticality conditions 

b 2 = 3.02-10‘3 cm-2 
prompt 

4. Radial extrapolation distance 
Sj. = (7.0±0.3) cm. 

5. Axial extrapolation distance 

= (8.5±0.3) cm. 

6. Infinite media multiplication constant 

= 1.1410.01 

7. Mean prompt neutron life time 
T = (70±4)*10'^ sec. 

8. Migration area for radial direction 

m 2 =(41±2)cm2 
r 

9. Migration area for axial direction 

m 2 =(45±2) cm2 

n 

10. Migration area anisotropy ratio 
A = l.lOlO.Ol 

11. Effective fraction of delayed neutrons 
iQeff = 0.0072+0.0003 

12. Temperature coefficient of reactivity 
|fl220C = l-5'10-4 (degree)-l 

The pulse-source method has some advantages which secure its success in the study of 
uranium-water systems. 

1. Measurement of the relation of cz(B^) within a rather wide range of the variation of 
and a very high accuracy of experimental value of a permit to determine a complex of the 
parameters M^, T with reliability and high accuracy. 

2. The assurance of the results produced is increased due to the absence of substantial 
corrections. The only significant correction on the variability of the process is important in 
the calculation of the highest moments of f (B^) — function. 

3. The precision of the results is determined particularly by the precision of the effective 
fraction /Sgf of delayed neutrons and by the precision of the effective reflector addition 
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This method provides with the principal possibility of the determination of S. 

4. The pulse-source method provides the more sensitive way of determination of the coef- 

ficient of anisotropy for the migration area A =_ji and of the dependence of the effective 

m2 

r 

addition on curve-radius of the boundary. 

5. The measurements with the calibrated instruments are simple and safe, they take rather 
little time since one deals with highly subcritical assemblies. These factors are substantial 
particularly in the examination of the series of lattices. 

6. The method permits to define the highest function-moments f (b 2). For this purpose it 
is essential to eliminate the correction effect on the variability and to reduce the effect of 8 
on the value of the geometrical parameter. It is practicable for the bare homogeneous systems. 

7. The pulse-source method can be used in the estimation of the critical dimensions for 
the instantaneous neutrons by the measurement of a = f (B^) for the assemblies with the critical 
mass of 0. 1-0.3, with the extrapolation of the dependence found of 

a = f (b 2) to the point a = 0. 
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FIG.A-1. TYPICAL CELL CONFIGURATIONS. 1 . (37-7) GAPLESS CELL. (61-1) GAPLESS 
CELL. 3 . (37+II-7) CELL. 4. (37-7) CELL WITH SINGLE ROW GAP 


5^,10^ cm~^ 




FIG.A-2. MATERIAL BUCKLING vs. BORON 
ROD BLACKNESS. STEEL CLAD FUEL ELE- 
MENTS, 8 MM SPACING, (37-7-t-n) CELL, 

1. CALCULATED RESULTS FOR 14xl2MM 
DIA. RODS. 2. CALCULATED RESULTS FOR 
10x8 MM DIA. RODS •, 0 - EXPERIMENTAL 
VALUES BASED ON TWO-GROUP CRITICAL 
EQUATION FOR TWO-ZONE REACTOR WITH 
RODS OF 14x12 MM DIA. AND 10x8 MM DIA., 
CORRESPONDINGLY. ■, 0 - AS ABOVE, 
BASED ON ONE GROUP APPROXIMATION. 


FIG.B-1. PARAMETERS OF THE LATTICES 
STUDIED vs. WATER TO URANIUM DIOXIDE 
VOLUME RATIO VH 20 /VUO 2 


- EXPERIMENTAL CURVE; 

CURVE 


- CALCULATED 


35e 


— 17 — 


Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 001 10006-3 





Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 0011 0006-3 



FIGXI-l . HIGHER HARMONIC CONTRIBUTION 
FOR VARIOUS SOURCE AND DETECTOR PO- 
SITIONS (r^.z, ARE SOURCE COORDINATES; 
tj, zj ARE DETECTOR COORDINATES). 

CURVE 1. r, = 0; Z, = 0.5H; rj = 0.54 R; 

Zd = |h. curve 2. r, = 0; Z, = 0.5 H; tj=0.54R; 
Zd = 0.8H. 

CURVE 3. Is = 0.44R; Z, = 0.5H; = 0.54R; 








FlG.C-2. THERMAL NEUTRON DENSITT' DECAY CONSTANT vs. 
BUCELING. CURVE 1. CADMIUM SHIELDED .ASSEMBLIES, Eq.(C--'O) 
CURVE 2. ASSEMBLIES SURROUNDED BY VATER REFLECTOR 

Eq.(C-I9). 



Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 001 10006-3 





